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In this paper we propose indirect probes of the scale of supersymmetry breaking, through obser- 
vations in the Extreme Universe Space Observatory onboard Japanese Experiment Module (JEM- 
EUSO). We consider scenarios where the lightest supersymmetric particle is the gravitino, and the 
next to lightest (NLSP) is a long lived slepton. We demonstrate that JEM-EUSO will be able 
to probe models where the NLSP decays, therefore probing supersymmetric breaking scales below 
5 x 10 6 GeV. The observatory field of view will be large enough to detect a few tens of events per 
year, depending on its energy threshold. This is complementary to a previous proposal where it 
was shown that 1 Km 3 neutrino telescopes can directly probe this scale. NLSPs will be produced 
by the interaction of high energy neutrinos in the Earth. Here we investigate scenarios where they 
subsequently decay, either in the atmosphere after escaping the Earth or right before leaving the 
Earth, producing taus. These can be detected by JEM-EUSO and have two distinctive signatures: 
one, they are produced in the Earth and go upwards in the atmosphere, which allows discrimina- 
tion from atmospheric taus and, second, as NLSPs are always produced in pairs, coincident taus 
will be a strong signature for these events. Assuming that the neutrino flux is equivalent to the 
Waxman-Bahcall limit, we determine the rate of taus from NLSP decays reaching JEM-EUSO's 
field of view. 

PACS numbers: 14.80.Ly, 12.60. Jv,95.30.Cq 



I. INTRODUCTION 

Probes of the origin and stability of the weak scale 
and proposed solutions to the standard model hierarchy 
problem is now underway at the Large Hadron Collider 
(LHC). On the theoretical side, the supersymmetric the- 
ory (Susy) arises as a solution with no significant devi- 
ation from the standard model (SM) in relation to elec- 
troweak precision observables. It has been shown [TJ [2] 
that neutrino telescopes can probe the scale of super- 
symmetry breaking (yF) or of Universal Extra Dimen- 
sions scenarios 3 . While R parity symmetry ensures 
that the lightest supersymmetric particle (LSP) is sta- 
ble, the identity of the LSP is determined by y/F. If 
Susy is broken at \/F > 10 10 GeV the LSP is typically 
the neutralino, if V F < 10 10 GeV it is typically the grav- 
itino. 

There are many Susy scenarios where y/F is low, 
among which Gauge Mediation Susy Breaking models [4] . 
In these scenarios the NLSP is a charged slepton, typi- 
cally a right-handed stau, and it's lifetime is given by pQ: 

-= ( -£L V ( 100GeV V 10 km , (1) 
^lO 7 GeVJ \ m fR J w 

where rrif R is the stau mass. It was shown pQCI] that Km 3 
neutrino telescopes can directly probe the Susy break- 
ing scale, when one consider scenarios where 5 x 10 6 < 
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\fF < 10 s GeV. In these scenarios, NLSPs produced in 
very high energy collisions will travel very long distances 
before decaying. A detectable flux of NLSPs can be pro- 
duced by the interaction of the diffuse flux of high energy 
neutrinos with the Earth. 

Here we consider [5] a complementary Susy breaking 
scale region, with y/~F % 10 7 GeV, which implies that 
the NLSP will decay after a short travel. After being 
produced by neutrino interactions in the Earth, a good 
fraction of these particles will decay inside the Earth or in 
the atmosphere, after escaping the Earth. In this work we 
consider NLSPs decay in the atmosphere (or right before 
reaching it). In a complementary investigation [51 [5] we 
analyze the y/~F region where they decay inside the Earth 
and might be detected by multi-km 3 neutrino telescopes. 
Once the NLSPs decay, taus (rs) will be produced and 
can be detected by fluorescence telescopes. We show that 
the JEM-EUSO observatory [7] can probe these scenar- 
ios, where NLSP decays will yield a few tens of events 
per year, depending on its lower energy threshold. 

We do consider the regeneration process, where ts de- 
cay into r neutrinos which in turn will charge current in- 
teract inside the Earth producing new ts. However these 
suffer large energy degradation, and only a small fraction 
of these events reach the detector with significant energy. 

As in PP, the crucial idea in probing \fF is that 
although the NLSP production cross section is much 
smaller than the one for SM lepton production, the NLSP 
range in the Earth is much larger than for a standard 
lepton. As shown in [2J, the NLSP energy loss is much 
smaller than the one for muons or ts. Here, when NLSP 
decays are considered, the NLSPs can be produced even 
further away from the detector, which implies that an 
extra decaying volume is gained. Also, as the NLSPs 
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are always produced in pairs, a considerable fraction of 
the observable events will consist of a pair of rs in the 
atmosphere. These coincident rs will be a distinctive 
indirect signature of NLSPs. Another strong signature 
comes from the fact that these ts emerge from the Earth 
and go upwards in the atmosphere, which discriminates 
them from down going atmospheric rs. 

In this paper we describe our analysis, where we de- 
veloped a Monte Carlo simulation to investigate r events 
generated from NLSP decays. The first steps, reviewed 
in the next section, reproduce the analysis done in [2], 
where NLSP production, propagation and energy loss are 
described in details. Subsequently, in section [TTTJ we de- 
scribe our simulation of NLSP decays and produced r 
propagation. The signatures and rates of these events 
in the JEM-EUSO observatory are determined in see- 
as well as their discrimination from the back- 
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ground. Finally we discuss our results and state conclu- 
sions. 



II. NLSP PRODUCTION 

Here we consider NLSP production by a diffuse flux of 
high energy neutrinos interacting in the Earth. The neu- 
trino flux as well as the NLSP production cross section 
and propagation are determined as described in detail in 
[2], and are used as the first steps in our Monte Carlo 
simulation. In the next section we describe the original 
part of our work, where NLSP decays are included, and 
the produced r rates in the JEM-EUSO observatory are 
determined. 

Although the diffuse flux of high energy neutrinos that 
reaches the Earth is yet unknown, there are several es- 
timates of its upper limit. Waxman and Bahcall (WB) 
[8] determined such a limit based on the observed cosmic 
ray flux, since neutrinos are produced from pion decays, 
which in its turn are produced from proton interactions. 
Considering optically "thin" sources, which would allow 
most of the protons to escape, they determine the neu- 
trino spectrum as 
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where the allowed interval depends on the cosmological 
evolution of the sources. Here we adopt the upper value 
of the WB limit as the cosmological neutrino flux that 
reaches the Earth. Our results can be translated to other 
neutrino fluxes by properly rescaling the WB limit. The 
initial neutrino flux contains both muon and electron 
neutrinos, in a 2 : 1 ratio. However we note that the 
initial neutrino flavor does not alter our results, since its 
interaction will always produce a left-handed slepton that 
will always immediately decay, having the right-handed 
slepton as a final product. For the same reason a large 
mixing of the cosmogenic neutrino flux does not modify 
our results. 



Once the neutrino flux is defined, the NLSP produc- 
tion cross section should be determined. We follow the 
cross section calculation described in detail in [2] and re- 
produce their results. In short the vN — > l^q process is 
analogous to the SM charged current interaction with q 
being a up or down type squark and, in the t-channel, 
the mediator is the chargino. The sub-dominant process 
with a neutralino exchange is also included in the cross 
section calculation. As a result of this interaction a II 
and a q will be produced. These will immediately de- 
cay in a chain that will always end with the production 
of two NLSPs, typically the right-handed stau (tr). We 
give our results considering the mass of the chargino and 
of the left handed slepton respectively as n%w = 250 GeV 
and m~ lL = 250 GeV, of the NLSP as m~ TR = 150 GeV 
and three possibilities for the squark mass, rriq = 300, 
600 or 900 GeV. 

The probability of a neutrino interacting in the Earth 
as well as the NLSP propagation through this medium, 
depends on the Earth density profile model. We use the 
model 9: described in detail in [TDJ. The lepton produc- 
tion cross section, from neutrino nucleon interaction is 
also determined in [ID] . 

Once the NLSP is produced it will propagate through 
the Earth and lose energy. The main process for NLSP 
energy degradation is photo-nuclear interactions. How- 
ever, as shown in [2] [IT], all radiative losses are sup- 
pressed due to the NLSP heavier mass when compared 
to standard leptons. 

NLSPs will always be produced in pairs and therefore, 
if one assumes the y/~F region where they do not decay, 
they will have a very distinctive signature in neutrino 
telescopes. As a cross check of our simulation, we re- 
produced the NLSP rate and energy distribution in Km 3 
neutrino telescopes as determined in [TJ [5] . 



III. 



NLSP DECAY 



The NLSP survival probability 

P s (x) — exp(—mx/Efcr) 



(3) 



is shown in Figure [I] as a function of the traveled distance 
x, for different neutrino energies E v and y/~F = 10 6 GeV, 
where the initial NLSP energy is typically E v /6 0. The 
NLSP production energy threshold is ~ 10 6 GeV, and 
defined by the q and the left-handed II- At these ener- 
gies, the distance jct is determined by Equation [I] and 
NLSPs decay when VF < 5 x 10 6 GeV. This feature is 
shown in Figure [5j The decay channel is tr —¥ r + G, 
where G is the gravitino and its mass is much smaller 
than the r mass. 

We developed a Monte Carlo simulation, where i9(10 5 ) 
events where generated, corresponding to neutrinos 
reaching the Earth, isotropically distributed both in en- 
ergy and impinging direction. These events are normal- 
ized by the WB limit. The NLSP production, propa- 
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FIG. 1. NLSP survival probability as a function of travel dis- 
tance, for different neutrino energies (as labeled) and y/F = 
10 6 GeV. 
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FIG. 2. Fraction of NLSP events in JEM-EUSO's field of view 
(FOV) versus y/~F for m,q reference values. NLSPs decay for 
VF < 5 x 10 6 GeV. 



gation and energy loss were simulated, reproducing the 
analysis shown in [2] (see Section [TTj) . 

The NLSP decay point is chosen from the decay proba- 
bility distribution (essentially 1 — P s ), and the generated 
r center of mass isotropic angular distribution is boosted 
into the laboratory frame. To a good approximation, it 
will follow the same direction as its parent NLSP. As two 
NLSPs are always generated in pairs, two rs will be pro- 
duced. A fraction of these ts subsequently decay, always 
generating a v T , and in 18% of these decays a muon is also 
produced. Once a v r is produced it can charge current 
interact producing a new r. We determine the fraction 
of the t energy carried by the v T as in [12, 13J. This 
regeneration process can happen a few times depending 
only on the r energy degradation. As mentioned before, 
this process will degrade the r energy and most of them 



FIG. 3. Representation of JEM-EUSO detection volume as 
a frustum of a cone, with a 100 Km height corresponding to 
the atmosphere altitude. NLSPs originated from neutrino in- 
teractions will decay in the atmosphere, producing observable 
rs. 

will not be detectable. 



IV. SIGNATURES AND EVENT RATES IN 
THE JEM-EUSO OBSERVATORY 

In order to determine the feasibility of NLSP indirect 
detection, where ts originated in NLSP decays would be 
probed by a large fluorescence telescope, we have to con- 
sider basic detector features. The JEM-EUSO telescope 
17] will observe fluorescence light emitted from an exten- 
sive air shower created by the interaction of a high energy 
particle in the atmosphere. It will orbit the Earth at an 
altitude of about 430 Km, yielding a detection area of 
~ 2 x 10 5 Km 2 with a 250 Km circular radius at the 
Earth surface. It is scheduled to be launched in 2016. 

At these energies, NLSP produced ts can be observed 
by JEM-EUSO mainly through the shower created once 
they decay. NLSP direct detection is hard, since the 
energy of its emitted fluorescence light is much less than 
the detector threshold, which is projected to be around 
10 19 eV. 

In order to determine the rate of observable events, 
we approximate the JEM-EUSO detection volume as a 
frustum of a cone, which represents the detector held of 
view. This is shown in figure [3] The frustum's height 
corresponds to the atmosphere altitude which is taken as 
100 Km. 

As described in the previous section, we simulate NLSP 
production from neutrino interactions in the Earth, their 
propagation and energy loss and finally their decay. The 
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FIG. 4. NLSP coincident decays in JEM-EUSO FOV as a 
function of \AF. 



neutrino energy ranges from 



10 b to 



where the lower limit corresponds to the NLSP produc- 
tion threshold. As the NLSPs are always produced in 
pairs, two rs will be produced from their decay. Al- 
though the pair of NLSPs travel in parallel, they can 
decay at different times, and each r can be produced at 
different positions. Each production point is chosen from 
the decay probability distribution. 

The r energy loss considers both ionization and ra- 
diation processes, where the latter includes loss due to 
bremsstrahlung, pair production and photo-nuclear in- 
teractions [21 [14] . We also follow the NLSPs which do 
not decay inside the Earth and reach the atmosphere, 
computing their energy and flux. The NLSP energy loss 
to the atmosphere is negligible. 

In summary we compute the r production point and 
initial energy distribution, as well as its decay point and 
energy distribution. These yield the number of rs which 
decay in or propagate through the detector field of view, 
and corresponding energy distributions. We also com- 
pute the number of coincident r decays in the field of 
view. All results are normalized by the WB neutrino 
flux. 

Our results are presented as a function of \/~F. Fig- 
ure [4] shows the number of coincident NLSP decays per 
year in JEM-EUSO field of view (FOV), while Figure [5] 
shows the NLSP average decay energy. As can be seen, 
the number of events is maximized for \JF = 3 x 10 6 GeV 
for both 300 and 600 GeV squark masses and 2 x 10 6 GeV 
for rriq = 900 GeV, due to its higher production energy 



threshold. However at these values of y/F, the NLSP 
average decay energy is lower than the detection thresh- 
old, and increases for lower y/~F values. We will further 
discuss this issue. 

Figure [6] shows the number of events per year in JEM- 
EUSO FOV as a function of the parent neutrino energy, 
for Vf-5x 10 6 GeV which maximizes the number of 
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FIG. 5. Average energy of events shown in Figure [4] Error 
bars represent standard deviation of each energy distribution. 
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FIG. 6. Events per year in JEM-EUSO FOV. Here m§ = 
300 GeV and \fF = 3 x 10 6 GeV, such that the number 
of events is maximized. The circles correspond to NLSPs 
(represented as f ) and triangles to rs. The NLSP (r) that 
decay are in yellow (green), the ones that go through the FOV 
without decaying in light blue (black), the ones which both 
NLSPs produced as a pair (or both rs) decay in the FOV in 
blue (purple). As a reference the total number of NLSPs that 
are produced in the direction of the FOV are shown in red 
and the WB limit as a black line. Oth 
have similar curves but with less events. 



events. It shows the NLSPs (f ) and rs that decay, the 
ones that go through the FOV without decaying, the ones 
which both NLSPs produced as a pair or both rs decay 
coincidently in the FOV. As a reference, the figure also 
shows the total number of NLSPs that are produced in 
the direction of the FOV and the WB limit. This figure 
shows results assuming rriq = 300 GeV, where for the 
other rriq values the shape of the curves are very similar 
but present less events. The total rate of events in JEM- 
EUSO per year, for the three values of rriq, are shown in 
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FIG. 7. Total number of r decays in JEM-EUSO FOV, re- 
quiring a minimal energy of 10 s GeV (top plot); 10 9 GeV 
(middle plot) and 5 x 10 9 GeV (bottom plot). 



TABLE I. Number of NLSP events per year in JEM- 
EUSO FOV, where %/F = 3(2) x 10 6 GeV for m 9 - = 
300,600(900) GeV. For m,- = 300 GeV these events corre- 
spond to the ones shown in Figure[6] "Total" is the total num- 
ber of NLSPs produced in the direction of the FOV; "Straight 
Through" refers to the NLSPs that go through the FOV; "De- 
cay " to those which decay in the FOV; and "Pair Decay" to 
the ones where the NLSP pair decays in coincidence in the 
atmosphere." 



rriq Total 


Straight Throuj 


;h Decay 


Pair decay 


300 | 1 x 10 7 


7, 2 x 10 5 


7,4 x 10 5 


2, 2 x 10 s 


600 1 1,3 x 10 6 


2,0 X 10 5 


1,3 x 10 5 


3,0 x 10 4 


900 1 2, 8 x 10 5 


1,7 x 10 4 


2,4 x 10 4 


7, x 10 3 



Table D 

As seen in this table a huge amount of NLSPs should 
either go through or decay in JEM-EUSO FOV. As the 
NLSPs arc not detectable due to the low amount of flu- 
orescence emitted while they transverse the atmosphere, 
the hope of probing these events relies on the rs pro- 
duced in their decay. However, although more than 10 5 
rs will tranverse or decay in JEM-EUSO FOV, most of 
these events have energies below the detector threshold 
(~ 10 19 GeV) and will not be observed unless this thresh- 
old is lowered. 

In order to determine the event rate that can be ob- 
served with an energy threshold around the currently 
projected by JEM-EUSO, we selected events above ar- 
bitrary lower energy values of 10 8 , 10 9 and 5 x 10 9 GeV. 
These are shown in Figure [7] and the integrated number 
of events, maximized as a function of y/F, in Table [ill 

As can be seen from Table [IT] there is still a reason- 
able amount of events for a maximized y/~F around JEM- 
EUSO energy threshold. Above 10 18 eV, 49 (19, 9) 
rs can be detected per year, for a 300 (600, 900) GeV 
squark mass. From these 5 (2, 1) r pairs decay in coinci- 
dence in the atmosphere in coincidence, which provides 
an excellent discriminating signature. For a threshold of 
5 x 10 18 GeV, a bit below the nominal projected thresh- 
old, 1 (0.4, 0.2) events per year can be observed by JEM- 
EUSO. 



A. Backgrounds 

Under the scenarios we are considering, NLSP decays 
will have a very distinctive signature. Since the NLSPs 
are produced inside the Earth and decay into rs that go 
upwards in the atmosphere, they can be discriminated 
from the more abundant down going cascades produced 
by normal ultra high energy cosmic rays. 

Also a considerable fraction of the NLSP pairs decay 
in coincidence in the atmosphere. Figure [8] shows the 
distribution of the time delay between two r decays in 
the atmosphere. These are shown for different \fF and 
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TABLE II. Number of r events per year in JEM-EUSO FOV 
after requiring a minimum r energy E T . %/F is such that it 
maximizes the number of events. The line corresponding to 
s/F — 3 x 10 6 GeV is shown for reference and has no energy 
requirement. 
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FIG. 9. Distance between two r decays in the atmosphere, 
assuming same parameters as in the previous figure. 



coincident showers can be well distinguished. 
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FIG. 8. Time delay between two r coincident decays in the 
atmosphere. Shown for arbitrary yF and lower NLSP energy 
requirements. The squark mass is 300 GeV but it does not 
vary considerably for larger values. 



energy lower limit requirements. It shows results for a 
300 GeV squark mass, which do not vary significantly for 
larger values. The time delay spread as a function of yF 
is small, and is distributed around 10~ 4 s. Considering 
that an atmospheric shower takes about 0.3 seconds, and 
that the detector resolution is about 2.5 /is, it is possible 
to determine the coincidence of these decays. 

Figure [9] shows the distance between the coincidence 
showers, considering the same parameters as for the pre- 
vious figure. Given that the detector spacial resolution 
is about 0.75 Km, it is clear that the majority of the 



V. DISCUSSION AND CONCLUSIONS 

In order to understand the effect of the experimen- 
tal energy threshold, we showed our results for arbitrary 
lower energy cuts. As shown in Table [TTJ , the lower 
threshold has a significant effect on how well large fluores- 
cence telescopes can probe the supersymmetry breaking 
scale. 

While an experimental energy threshold of 10 18 eV al- 
lows JEM-EUSO to discover NLSP produced rs and con- 
sequently probe a significant y/F region, a 10 19 eV allows 
to set limits in the production and \/F parameters. Ta- 
ble n] shows that, for a maximum value of VF, i?(10 5 ) 
of detectable events will go through JEM-EUSO FOV. 
From these hundreds of events per year can be seen if 
the energy threshold is set to 10 17 eV, while tens for 

10 18 and more than a year is needed to detect events for 

10 19 eV. Although these numbers are for a maximized 
value of yF, one can see from Figure M that there is a 
significant range, in which 10 5 j£ VF % 10 6 GeV, that 
can be probed with the same order of events. 

We have shown that large fluorescence telescopes such 
as JEM-EUSO have the potential to indirectly detect NL- 
SPs which are modeled in various supersymmetry break- 
ing scenarios. We show that scenarios where the su- 
persymmetry breaking scale vF is such that the NLSP 
will decay close to or in the atmosphere can be probed 
by JEM-EUSO. Depending on the experimental en- 
ergy threshold, scenarios with 10 5 ^ \f~F % 10 6 can be 
probed, and NLSPs can be indirectly searched. 

This work complements the direct probe for long-lived 
NLSPs [HE], where scenarios with 5 x 10 6 < VF < 5 x 
10 8 can be directly probed by neutrino telescopes. It is 
also complementary to searches for NLSP decays at the 
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